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Studies of Majorana bound states in semiconducting nanowires frequently neglect the orbital effect 
of a magnetic field. Systematically studying its role leads us to several conclusions for designing 
Majoranas in this system. Specifically, we show that for experimentally relevant parameter values the 
orbital effect of a magnetic field has a stronger impact on the dispersion relation than the Zeeman 
effect. While Majoranas do not require the presence of only one dispersion subband, we observe 
that the size of the Majoranas becomes unpractically large, and the band gap unpractically small, 
when more than one subband is filled. Since the orbital effect of a magnetic field breaks several 
symmetries of the Hamiltonian, it leads to the appearance of large regions in parameter space with 
no band gap whenever the magnetic field is not aligned with the wire axis. The reflection symmetry 
of the Hamiltonian with respect to the plane perpendicular to the wire axis guarantees that the wire 
stays gapped in the topologically nontrivial region as long as the field is aligned with the wire. 

PACS numbers: 73.63.Nm, 74.45.+C, 74.78.Na 


I. INTRODUCTION 

The search for Majorana bound states, the simplest 
non-Abelian particles, is fueled by their suitability for 
fault-tolerant quantum computation.^’^ A large fraction 
of the experimental effort^“^ is focused on creating Ma¬ 
joranas in semiconducting nanowires with proximity su¬ 
perconductivity, spin-orbit coupling, and magnetic field. 
The theoretical foundation for this platform was initially 
developed for a single one-dimensional spinful band with 
intrinsic superconducting pairing.Due to its compact¬ 
ness this model can be solved analytically, and it predicts 
that Majorana bound states appear when -f , 

when the Zeeman energy becomes larger than the har¬ 
monic mean of the superconducting gap and the chemical 
potential. 

The single-mode model is minimalistic and neglects 
many physical phenomena that are crucial for under¬ 
standing the properties of the Majorana bound states. 
The existing extensions of this model study multimode 
wires,better modeling of the induced gap,^^’^^ the role 
of electrostatics,^^ disorder,and the k -p-model.^^ 
The orbital effect of a magnetic field was analyzed both 
in planar wires^^’^^ and on the surface of a cylinder. 

We systematically study the influence of the orbital 
effect of a magnetic field on the symmetries of the Hamil¬ 
tonian and the topological phase diagram for a three- 
dimensional (3D) nanowire. The orbital effect of a mag¬ 
netic field perpendicular to the wire induces a skipping or¬ 
bit motion of the electrons. The cyclotron radius becomes 
comparable to the typical wire diameters d ^ 100 nm 
already at the field of 0.3 T, and at chemical potential 
corresponding to the optimal topological band gap. In 
addition, a field parallel to the wire shifts the energies of 
each band due to the effect of magnetic flux. We expect 
the shift of the energies to be comparable to the level 
spacing when the flux through the wire diameter is of the 
order of a flux quantum. Our findings are very different 


from those of Refs. 18-20 because we do not limit our 
analysis to a Hamiltonian with an artificially high spatial 
symmetry, or low dimensionality. 


II. MODEL 

We consider a 3D semiconducting nanowire with 
Rashba spin-orbit coupling and proximity-induced s-wave 
superconductivity. The nanowire cross section is a regular 
hexagon, and the nanowire is translationally invariant in 
the x-direction. Its Bogoliubov-de Gennes Hamiltonian 
(BdG) is 

-ffBdG = - lfjT^+a [pyGx - PxCTy) Tx 

+ +ATj,, (1) 

and it acts on the spinor wave function At = 

, where 'tpe, V’h are its electron 
and hole components, and are the spin-up and 

spin-down components. We introduced the Pauli ma¬ 
trices ai acting on the spin degree of freedom and 
acting on the electron-hole degree of freedom. Fur¬ 
ther p = —ihV + eAxz is the canonical momentum, 
with e the electron charge and the vector potential 
A = [By{z - Zo) - B^{y- yo), 0, B^iy - yo)f chosen 
such that it does not depend on x. We set the offsets yo 
and zq to ensure that the average vector potential vanishes 
in the superconductor. This choice corresponds to a limit 
when the superconductor is thinner than the screening 
length and its total supercurrent is zero, appropriate for 
existing devices. Finally, m* is the effective electron mass, 
Ez = /2 the Zeeman energy, A the superconduct¬ 

ing pairing potential, a the Rashba spin-orbit coupling 
strength, and y = yo + the chemical potential created 
by a constant electric field E in the sample parallel to 
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the z-axis, such that the Rashba spin-orbit acts in the 
x^-plane. 

First we consider a model with a constant superconduct¬ 
ing gap A inside the wire [see Fig. 1(a)] and then proceed 
to make a more realistic model of the superconductor. 
To do that we set the superconducting order parameter 
A to zero in the wire and add a superconductor to the 
top which covers 3/8 of the circumference of the wire [see 
Fig. 1(b)]. We choose the thickness of the superconductor 
to be 20 nm and set A in the superconductor such that 
the induced gap of the lowest band is Aind = 0.250 meV. 
This is done by computing band energies at /c = 0 over 
a range of /i and matching the minimum to Aind- We 
add a tunnel barrier between the two materials to change 
the transparency of the superconductor. In the setup of 
Fig. 1(c), we break the reflection symmetry with respect 
to the x 2 ;-plane by moving the superconductor to the side 
similar to the experimental setup of Mourik et al.^ 

To perform the numerical simulations we discretize 
the Hamiltonian on a cubic lattice with lattice constant 
a = 10 nm, much smaller than the minimal Fermi wave¬ 
length in the parameter range we consider. The dis¬ 
cretization does not break or introduce any additional 
symmetries. The Hamiltonian at a lattice momentum 
k equals H {k) = h texp(i/c) -h exp(—i/c) where h is 
the Hamiltonian of the cross section of the tight-binding 
system and t is the hopping matrix between the neighbor¬ 
ing cross sections. We introduce the vector potential by 
Peierls substitution tnm tnm exp(—ie j Adl)?^ We per¬ 
form the numerical simulations using the Kwant code.^^ 
The source code and the specific parameter values are 
available in the Supplemental Material.The resulting 
raw data are available in Ref. 24. 


III. SYMMETRY ANALYSIS 

The Majorana bound states are protected by the com¬ 
bination of the band gap and the particle-hole symmetry 
of the Hamiltonian VH {k) V~^ = —H {—k). In the ba¬ 
sis of Eq. (1) this symmetry has the form V = cr^Ty/C, 
with 1C the complex conjugation. In general there are no 
additional symmetries and the Hamiltonian belongs to 
symmetry class Particle-hole symmetry only requires 

that the energy En{k) of n-th band at momentum k is 
En{k) = —Em{—k) of some other m-th band; at the same 
time V puts no constraints on E^ itself. This means 
that whenever E^ changes sign at a certain momentum, 
the band structure becomes gapless. This tilting of the 
band structure^^ [shown in the middle panels in Fig. 2, 
where En{k) ^ —Em{k)] is a strong effect that does not 
vanish with superconducting gap or spin-orbit, and can 
easily become larger than the induced gap, rendering the 
creation of Majoranas impossible. 

The tilting of the band structure is absent if the Hamil¬ 
tonian has an extra chiral symmetry alongside V. It has 
been shown that the Hamiltonian has an approximate 
chiral symmetry CH (k)C~^ = —H (/c), C = valid 




Figure 1. Three hexagonal nanowire devices we consider: (a) 
(a) with an intrinsic pairing term and with a proximity-coupled 
superconductor (b) on the top and (c) on the side. The last 
two setups have tunnel barriers between the superconductor 
and the nanowire. 


when the wire diameter d is smaller than the spin-orbit 
length Iso = k?/nnCa and By = 0. Then the PyCx^z 
term, associated with the transverse motion in Eq. (1), is 
negligible. Without the tilting, the system is gapped in 
every region of parameter space, except at the topological 
phase boundaries. However, for relevant experimental 
parameters,^ the orbital terms break this symmetry more 
strongly than the spin-orbit term, bringing the system 
back to symmetry class D. 

We perform a systematic search of symmetries that 
the Hamiltonian (1) may have.^^ We find the reflection 
symmetry with respect to the ^ 2 :-plane IZxH (k) 1Z~^ = 
H (— /c), IZx = cFx^{x + x'). It is independent of the 
wire geometry and spin-orbit strength and guarantees 
the absence of tilting whenever the field is aligned with 
the x-axis. The combined symmetry V' = ExV is local 
in momentum space and ensures the absence of band 
structure tilting: V'H{k)V'~^ = —H{k). 

Additionally, we find a chiral symmetry C' = TylZy^ 
C'H (k) = —H (/c), with IZy = (JyS{y + y') the reflec¬ 

tion with respect to the ^-axis. This chiral symmetry 
holds when the magnetic field lies in the xz-plane and 
none of the potentials in Eq. (1) break IZy, like in the 
setups of Eigs. 1(a) andl(b). When present, C' guar¬ 
antees the absence of band structure tilting just like C. 
This symmetry is present in most theoretical models, and 
in particular it is obeyed by the Hamiltonians used in 
Refs. 18-20. A finite By breaks both IZx and C' therefore, 
the bands can tilt and close the topological gap. The band 
structures in Fig. 2 summarize the relation between the 
geometry of the setup of Figs 1(b) and Fig. 1(c), magnetic 
field orientation, and the symmetries of the Hamiltonian. 


IV. CALCULATING THE TOPOLOGICAL 
PHASE DIAGRAM 

We use an optimized algorithm to quickly find all the y 
values corresponding to the topological phase transitions 
at once. The topological transitions in symmetry class 
D occur when Pf i^BdG(^ = 0) changes sign.^^ Since the 
sign change of Pf HbcIG is accompanied by the appearance 
of zero energy states, we need to find y and such that 
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Figure 2. Band structures of the setup of Fig. 1(b) (top) and 
Fig. 1(c) (bottom). Each panel is labeled with the symmetries 
respected by the corresponding Hamiltonian. The dashed black 
line indicates the Fermi energy (E — 0). The red dashed lines 
show the size of the band gap if it is present. In the top row, 
the reflection symmetry of the wire along the y-axis IZy makes 
the Hamiltonian have a chiral symmetry C' when the magnetic 
field lies in the xz-plane. The wire used for the calculation 
of the bottom row dispersions lacks IZy and therefore has no 
C'. Without C' the bands are allowed to tilt and the gap may 
close whenever By ^ 0 or Bz ^ 0. A magnetic field parallel 
to the x-axis preserves IZx, which protects the band gap from 
closing. 

= 0 ) 7 /^ = 0. Using that /i enters the Hamiltonian 
as a prefactor of a linear operator, we rewrite this equation 
as a generalized eigenproblem: 

HBdG {fi = 0,k = 0)'tp = liT^tp. (2) 

The real eigenvalues of this eigenproblem are the values 
of ji where the gap closes at /c = 0 [see Fig. 3(a)], and 
they can be found using standard generalized eigensolvers. 
If the dispersion relation is gapped also at any finite /c, 
these gap closings are the boundaries of the topological 
phase. 

Since the eigenvalues of i^BdC come in opposite sign 
pairs, the real eigenvalues of Eq. (2) always come in de¬ 
generate pairs, and each pair lies at a transition between 
a trivial and a nontrivial phase. We complete the cal¬ 
culation of the topological phase diagram by using as 
a reference point that i^BdcC/^ = — 00 ) is topologically 
trivial. 

The generalized eigenvalue algorithm for finding phase 
boundaries does not guarantee that H{k) is gapped for 
/c 7 ^ 0 , and therefore we calculate the magnitude of the 
gap Egap in the topologically nontrivial regime separately 
for each set of parameter values. We form a translation 
eigenvalue problem to calculate all the modes of i^BdC 
at a given energy E and check whether there are any 
propagating modes.By using a binary search in E 
for the energy at which the propagating modes start to 
appear, we find Egap [see Fig. 3(b)]. 
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Figure 3. (a) Energy spectrum at /c = 0 of the setup of Fig. 1(a) 
as a function of chemical potential y. The blue points are the 
solutions of HBdG^ = E'lp at fixed y marked by the blue line. 
The green points are the real eigenvalues of Eq. (2) lying at 
E = 0 (the green line), (b) The gap size for the same setup 
and parameters, with dark gray regions trivial and the orange 
regions topological. 

The real space size of the Majoranas ^ imposes a lower 
bound on the nanowire length required to create them. 
To calculate ^ we find the eigenvalue decomposition of the 
translation operator at zero energy. The eigenvalue Amin 
closest to the unit circle corresponds to the slowest decay¬ 
ing part of the Majorana wave function. We calculate ^ 
using 

^ = |l0g“tAmin||- (3) 


V. RESULTS 

We use realistic parameters of an InSb nanowire a = 
20meVnm, m* = O.Olbmg, A = 0.250meV, d = 100nm, 
and ^ = 50. At the high fields that are typically used 
in experiment (E > 1 T), we find that the Zeeman effect 
of the magnetic field has a lower impact on the phase 
boundaries than the orbital effect of the magnetic field 
(see Fig. 4). We verify that the band gap is protected by 
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Figure 4. Phase diagrams of the setup of Fig. 1(a) (a), (b) 
without the orbital effect of a magnetic field and (c), (d) with 
it. The green lines depict the topological phase transitions. 
The colored regions are topologically nontrivial, with the color 
representing the size of the topological band gap F^gap- At 
B > IT the orbital effect of the magnetic held becomes 
stronger than the Zeeman effect and changes the sign of the 
slope of half of the phase boundaries. Furthermore, the orbital 
effect leads to a faster suppression of the band gaps with 
magnetic held. The narrow regions with suppressed F^gap 
originating from the crossings of the phase boundaries in (c) 
are due to Dirac cones appearing in (fc, S)-space and are 
protected by C' . The vertical black line in (a) indicates the 
value of the magnetic held used in Fig. 3. 



Figure 5. Same as Fig. 4, but with the magnetic held slightly 
misaligned. We observe that the band gaps close quickly 
upon changing the direction of the magnetic held towards the 
spin-orbit direction in y. 
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Figure 6. Phase diagrams of the setup of Fig. 1(c) (a) without 
the orbital effect of a magnetic held and (b) with it. Color 
scale corresponds to F^gap, with the topological regions colored 
and trivial regions in grayscale. The histograms in the right- 
hand panels show the distribution of the gap values sampled 
in the topological regime within the selected parameter range. 
Neglecting the orbital effect of the magnetic held incorrectly 
leads to a strong dependence of the critical held on fi. With 
the orbital effect of magnetic held hux penetration through 
the quasiparticle trajectory changes the interference phases 
and can suppress the topological gap F^gap- 


C' as long as = 0, despite that the orbital effect of the 
magnetic held reduces F^gap- 

In agreement with our expectations a hnite < 0.1 T 
leads to the closing of the band gap (see Figs. 2 and 
5). The maximum tolerable By becomes smaller with 
increasing ji. The narrow regions with suppressed F^gap 
visible in Figs. 4(c) and 4(d) are the consequence of Dirac 
cones appearing in {kx^ 5)-space and are protected by C'. 
Breaking IZy breaks C' and removes these Dirac cones. 

We now turn to study the system shown in Fig. 1(c) 
that has C strongly broken and only IZx and V remaining. 
Since the induced superconducting gap Aind ~ 250 peV 
in Ref. 3 is much smaller than the NbTiN gap 2meV, 
the system must be in the long junction limit, where 
F^Th ^ A. In the long junction limit the induced gap 
equals Aind ~ hTv-p/d^ where T is the transparency of 
the tunnel barrier, and v-p the Fermi velocity. In the ab- 
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the Majorana properties in the presence of an orbital field, 
we have repeated the calculations shown in Figs. 6 and 
7 using a fivefold larger spin-orbit strength reported in 
Ref. 31. We find that the topological band gap increases 
overall and in particular the maximal gap grows from 
0.14 meV to 0.21 meV, while the minimal decay length 
remains almost the same. Therefore, increasing spin-orbit 
strength has a positive but not very strong effect on the 
topological band gap. 


VI. DISCUSSION AND CONCLUSIONS 



Figure 7. Same as Fig. 6, but with color representing inverse 
Majorana length . The histogram and color scales are 
truncated from above at 1 pm“^. The mode of the distribution 
of reduces from 0.35 pm“^ to 0.10 pm“^ upon taking the 
orbital effect into account. Although the Majorana lengths 
are overall much larger with the orbital effect of the magnetic 
field, the minimal length is close to 200 nm in both cases. 


sence of the orbital effect of a magnetic field, this means 
that the Zeeman energy has to exceed Aind and therefore 
the critical value of the magnetic field at which the gap 
closes strongly depends on /r as seen in Fig. 6(a). With 
the orbital effect of the magnetic field fiux, penetration 
through the quasiparticle trajectory changes the interfer¬ 
ence phases, which suppresses the induced gap and causes 
the topological phase transitions to occur at a value of B 
corresponding to a single fiux quantum penetrating the 
wire area [see Fig. 6(b)]. 

The Majorana decay lengths ^ significantly increase 
when including the orbital effect of the magnetic field in 
the Hamiltonian (see Fig. 7). Specifically, the mode of 
the distribution of ^ changes by a factor of ^ 4 in the 
parameter range we consider (see histograms in Fig. 7). 
However, the minimum values of ^ without orbital effect 
and with it are both ^ 200 nm. Therefore, /i needs to be 
tuned with sub-meV precision within the lowest band in 
order to create Majorana bound states with practically 
relevant parameters. 

To investigate the effect of the spin-orbit coupling on 


We have shown that the orbital effect of a magnetic field 
complicates the creation of Majoranas in nanowires. Or¬ 
bital terms break the chiral symmetry C and prevent the 
appearance of Majoranas whenever the magnetic field is 
not aligned with the wire axis. When the field does point 
along the x-axis, we find that the refiection symmetry TZx 
in combination with particle-hole symmetry V protects 
the band gap from closing everywhere in (H,/i)-space, 
except at the topological phase boundaries. At exper¬ 
imentally relevant values of magnetic field, the orbital 
effect has a stronger impact on the dispersion relation 
than the Zeeman effect. Furthermore, the orbital effect 
suppresses F^gap and increases However, the maxi¬ 
mum value of the F^gap in the topologically nontrivial 
region does not change as drastically (from 0.21 meV to 
0.14 meV) and the minimal decay length changes even 
less (from 201 nm to 210nm). The refiection symmetry 
IZx of the Hamiltonian that we consider is respected by 
any Rashba spin-orbit interaction. Dresselhaus spin-orbit 
coupling breaks Tlx] however, it is expected to be weak 
in the nanowires. 

Our simulations can be made more complete by com¬ 
plementing them with self-consistent electrostatics and 
magnetic field screening by the superconductor. An addi¬ 
tional extension of our work is to go beyond the effective 
mass approximation and to use the k-p model. A separate 
topic of study is the interplay between the orbital effect 
of the magnetic field and disorder. We expect that the 
sensitivity to disorder will increase by taking the orbital 
effect of the magnetic field into account. 

Our results suggest that keeping the chemical potential 
low is required to obtain Majoranas with reasonable length 
and energy scales. Furthermore, our findings reveal a 
complication in realizing more sophisticated Majorana 
setups, such as a T-junction required for braiding. This 
is because of the requirement that the field should be 
aligned with the nanowire. A possible strategy to reduce 
the undesirable orbital effect of the magnetic field is to 
use nanowires with smaller diameters at a cost of reduced 
electric field effect and increased disorder sensitivity. 
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